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ABSTRACT 



Via joint analysis of a calibrated telescopic survey, which found scattering 
Kuiper Belt objects, and models of their expected orbital distribution, we mea- 
sure the form of the scattering object's size distribution. Ruling out a single 
power-law at greater than 99% confidence, we constrain the form of the size 
distribution and find that, surprisingly, our analysis favours a very sudden de- 
crease (a divot) in the number distribution as diameters decrease below 100 km, 
with the number of smaller objects then rising again as expected via collisional 
equilibrium. Extrapolating at this collisional equilibrium slope produced enough 
kilometer-scale scattering objects to supply the nearby Jupiter-Family comets. 
Our interpretation is that this divot feature is a preserved relic of the size dis- 
tribution made by planetesimal formation, now "frozen in" to portions of the 
Kuiper Belt sharing a "hot" orbital inclination distribution, explaining several 
puzzles in Kuiper Belt science. Additionally, we show that to match today's 
scattering-object inclination distribution, the supply source that was scattered 
outward must have already been vertically heated to of order 10°. 

Subject headings: comets: general — Kuiper belt: general 
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Introduction 



Measurements of the Kuiper Belt's size distribution (number at each diameter D) 
constrain accretional processes at planet formation and, potentially, subsequent collisional 
or physical evolution. Because astronomers observe brightnesses rather than D, object 
absolute magnitudes H are tabulated as the observable proxy for the size distribution. If 
the albedo p is known, H is convertible to D; for 5% albedo a D = 100-km object has 
if =9. 16 in g band. Using H g magnitudes avoids uncertainties introduced by assuming the 
unknown albedos. 

Collisional and accretional theories suggest exponential forms for the N(H) distribution. 
A differential number distribution of the form dN/dH oc lO"^ 9 with a logarithmic 'slope' 
a corresponds to a power- law D distribution dN/dD oc £)-( 5a + 1 ). Although power-law 
D distributions provide acceptable fits to Kuiper Belt surveys over spans of a few 

neces s ary over larger H ranges 



(Jewitt et al. 


1988; 


Gladman et al. 


2001 


Fraser & Kavelaars 


2008|). For the steep 



Bernstein et al. 



2004 



Fuentes &: Holman 



2008 



detections are dominated by objects near the largest H magnitude (smallest size) visible 
in a given survey. An a > 0.6 slope cannot continue as H — > oo (D — > km) or the total 
mass diverges; thus a slope change (generically called a break) is required. Evidence of 
such a break now exists for trans-Neptunian objects (TNOs) in the main Kuiper Belt (at 
distances d ~ 38-46 AU), both near the sensitivity limit of ground-based telescopic surveys 



( IFraser & Kavelaars 



from deeper HST (IBernstein et al. 



2008 



Fuentes fc Holman 



20081 ) (reaching H g ~9-10 at 40 AU) and 



20041 ) observations (H g ~13 at 40 AU). This break has 



been modelled as a gradual transition to a smaller value of a, a "knee". 

Probing beyond a break is difficult because small TNOs are faint. This problem is 
reduced when observing the scattering objects (SOs); these are mostly TNOs with perihelia 
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q < 35 AU (see below) and thus smaller objects are detectable as they come near the 
observer. At any time some SOs are only d = 20-30 AU away, allowing a 4-m telescope, in 
excellent conditions, to detect objects down to H g ~ 12. For a monotonically increasing 
number distribution N(H g ), the abundant small objects at the observable volume's 
innermost edge should dominate the detected sample. This is not what our survey found, 
necessitating a relative lack of small SOs. 



2. Models 

Several models of the SO orbital element distribution were exposed to the calibrated 
observational biases of the Canada France Ecliptic Plane Survey (CFEPS) in order to 
quantitatively constrain the intrinsic N(H g ) distribution. Drawing SOs from an orbital 



distribution mode 



survey simulator (IJones et al 



, and selecting H g from a candidate N(H g ) distribution, the CFEPS 



20061 ) determines each object's observability and produces a 



set of "simulated detections" expected from the model. 

Tw o different simulations of this population are from a modified version of 



KaibetaL 



( 1201 lbh (henceforth KRQ11). KRQ11 focuses on the effects of solar migration in the Milky 
Way on Oort Cloud structure. While this is not the focus of the current work, we can use 
the KRQ11 control simulations, which assume an unchanging local galactic environment 
like our current one. 



To test the sensitivity of our results to the choice of the dynamical contex t, we 



performed the same analysis on an independent model. 



Gladman fc Chan! (120061 ) modelled 



the scattering of objects an ini tial Solar System with an additional planet (of order Earth's 



size). As previously reported ( iPetit et al 



201 ll ) this model also (perhaps surprisingly) 



reasonably represents the current SO (a, q) distribution , although too "cold" in inclinations. 
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In fact, this model and the cold KRQ11 model produced very similar results. We find 
that the results for this model are nearly identical to KRQ11 results, showing that our 
conclusions are not strongly dependent on the model of Solar System history, as the objects 
that are currently in the Centaur and detectable SO region (mostly a <200 AU) today have, 
unsurprisingly, almost forgotten their initial state except for the the inclination distribution. 
Thus, the current SO orbital distribution is not diagnostic of the number and position of 
the planets early in the Solar System's history. 



CFEPS provi ded a set o 



calibrated survey (IJones et al 



Observations 



detections of outer Solar System objects in a precisely 



2006 



Kavelaars et al. 



20091 ) whose pointing history, detection 



efficiency, and tracking performance were recorded. The final set of TNO detections (with 
full high-precision orbits) and the fully calibrated pointing history make up the L7 release 



(IPetit et al. 



201f[ ). This absolute calibration of CFEPS allows a model of the present orbital 
(and size) distribution of to be passed through the CFEPS Survey Simulator, yielding a 
set of simulated detections whose orbital and H g distributions can be compared to the real 
detections. 

Our models provide orbital distributions of all TNOs. The "scattering" TNOs are then 
selected out of the final 10 Myr stage of the model integrations usin g the criteria: variatio n 



of a >1.5 AU in semim ajor axis during 10 Myr, with a < 1000 AU (IGladman et al. 



2008 



Morbidelli et al. 



( iDuncan fc Levison 



2004). H istorically, a simple q cut was used t o isolate the "scattered disk" 



1997 



Luu et al. 



1997 



Trujillo et al 



2000l ). which undesirably includes 



resonant objects and most inner main-belt TNOs. 



The CFEPS SO sample consists of 9 detected objects (Table 1). In addition, we 
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include two objects discovered in a high-latitude extension survey (covering ~ 470 sqdg in 
2007-2008, extending up to 65° ecliptic latitude), which was fully calibrated in the same 
way as CFEPS. 

The CFEPS Survey Simulator was used to observationally bias the three models. Given 
an orbital model and a candidate N(H g ), the simulator determines if one of the survey fields 
would have detected test model objects, producing a distribution of "detected" objects as if 
the model were the true orbital distribution. The resulting set of "observationally-biased" 
objects are statistically compared to the set of real detections. 

To characterize the form of the N(H g ), we introduce a novel parameterisation which 
allows for the exploration of distributions with knees and divots (a sudden drop in 
the differential number of objects followed by a recovery). We parameterised the size 
distribution (Fig. 2A) with the fixed slope = 0.8 (see below) for SOs brighter than a 
break at H g = 9 (D ~100 km), allowed an adjustable slope a/ for fainter objects, and an 
adjustable contrast c > 1. 

The if s -magnitudes are drawn from one of three types of distributions: 

1. a single exponential of logarithmic slope a, 2. a N(H g ) with a knee (contrast c =1). 
That is, one slope at> for SOs with H g < H knee and atf for H g > H knee , where the N(H g ) 
distribution is continuous across the knee at H krifr and negative slopes a/ are allowed as 



per suggestions in the literature ( IBernstein et al. 



20041 ). and 3. one slope to a divot 



at Hdivot, which is a sudden drop in differential number by a factor c, with a potentially 
different slope atf beyond the cliff at H = Hdmot- Although in reality the discontinuity is 
unlikely to be an instantaneous drop, our data do not merit trying to constrain the values 
of the expected steep negative slope and small ex tent over which it drops; collisional models 



( IFraser 



20091 



Campo Bagatin &: Benavidez 



20121 ) do show collisional divots where the drop 



occurs over short diameter range of much less than a factor of two in D (a few tenths of 
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magnitude in H g ). 



In principle there are four parameters: a&, a/, H divot , and c (Fig. 2). For H g < 9 a 
single power law of a& ~ 0.8 does match our detections, and we elected to fix this slope 
at that value with the unifying philosophy that all the hot transneptunian populations 
share this sa me hot slope; ah = 0.8 matches both the hot Classical belt, measured down 



to H g ~ 8.0 flPetit et al 



2011 



Fraser & Kavelaars 



2008 



Fuentes &: Holman 



20081 ). and 



20121 ). both of which 



to the 3:2 resonators, measured down to H g ~ 9.0 ( IGladman et al. 
appear consistent with single power laws for H brighter than those limits. Our detections 
require a transition around H g =9-10 to explain the lack of small detections; we thus fix 
the knee/divot for our analysis at H g = 9 (slightly larger than D=100 km for 5% g-band 
albedo). This leaves only two free parameters: the contrast c at the divot and the slope otf 
for objects with absolute magnitudes fainter than the divot/knee. 

To assess a match, the Anderson-Darling (AD) statistic is calculated between our 
11-object sample and the distribution of simulated detections from the model, for each 
orbital parameter. An AD significance level of < 5% means a rejection of the hypothesis 
that the real SO observations could be drawn from the simulated detections at the 95% 
confidence level (for that orbital parameter). For a model to not be rejectable we use the 
criteria that none of the q, d, i, and H g distributions can be rejectable at > 95% confidence. 



4. Absolute Magnitude Distribution 

The observational bias is strong (Fig. 1) , but when accurately calibrated allows us to 
constrain the H g distribution's form. Single power-laws predict significantly more close-in 
detections than were seen by CFEPS; for a slope of a = 0.8, roughly half of the expected 
detections (Fig. 1 D's blue dashed curve) should have a distance at detection d < 23 AU, 
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which is the closest real SO in our sample. This observationally biased model predicts that 
the majority of detected SOs would have orbits with q <20 AU at d — 20-25 AU and be 
small (H g > 9 or D < 100 km) objects, in contrast to our detections, which demonstrates 
that our observations are sensitive beyond the break. This discrepancy exists for all of our 
models when using a single power- law, and thus must arise from the assumed N(H g ). We 
rule out a single power-law at 95% confidence. 

Our relatively small sample is powerful because our detected SOs span the break, and 
when coupled with the precise CFEPS calibration, allows the non-detection of H g =10-12 
SOs (several magnitudes past the divot) to provide a strong constraint on N(H g ). CFEPS 
detected moving objects as close as 20 AU with no rate of motion dependence down to 
this limit. Because our orbits are accurate, we can separate the SOs from the other hot 
populations, and use a dynamical model specific to the SOs. 

In Fig. 1, our best N(H g ) case has the obvious and previously- known p roblem that 



the model's orbital inclinat ions are mostly lower than the true population's (IPetit et al. 



2011 



Gladman et al 



20121 ). Fig. 1 shows a divot (c ~ 6, af — 0.5) which produces a 



good match between the model's expect ed detections 



(red), excepting the % problem. Models (ILevison et al. 



green curve) and the real SO sample 



20081 ) which scatter out a cold TNO 



population, from d < 30 AU and with initial inclination distribution widths < 6°, to 
eventually form today's hot population produce current TNO populations where too many 
low-z detections are expected in observational surveys (Fig. IB). This is part of growing 
evidence that the original planetesimal disk supplying to day's high-i objects must have 



already been vertically excited be fore being scattered out, (IPetit et al. 



2012 



Brasser fc Morbidelli 



2011 



Gladman et al. 



20121 ) which ha s strong cosmogonic implications for an extended 



quiescent phase of the early Solar System ( jLevison et al 



20081 ). We therefore computed a 



new SO model with a hotter (cr, ~ 12°) initial disk; this provides an excellent match with 
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today's SO i-distribution (see Fig. 4) and we constrain N(H g ) below using this model. 



To constrain the size distribution, a grid of possible divot contrasts and post-divot 
slopes was explored. The parameter amongst i, q, d and H which is most rejectable is 
recorded; Fig. 3 shows acceptability levels for the range of explored parameter pairs (c, aj). 
A single power-law of a = 0.8 (blue star Fig. 3) has < 1% probability. We are left with a 
range of acceptable parameter space, including knee (c=l) and divot (c >1) scenarios; we 
further constrain N(H g ) by looking to other Kuiper Belt populations. 

The so-called hot Kuiper Belt populations (the hot main belt, i nner belt, resona nt, and 
detached TNOs) share an i distribution half-width of roughly 15° OPetit et al.l 120111 ) with 
the SOs, suggesting a cosmogonic link. In analyses of deep luminosity functions dominated 



2004 



Fuentes et al. 



by ho t main-belt detections, the common conclusion (IBernstein et al. 
20101 ) was that for magnitude g > 25 the slope must break to a faint a.f < 0.3 value or even 
become negative in order to explain the lack of detections in the following few magnitudes; 
beyond this no data exists for the main Kuiper Belt. For SOs and their companion objects 
(Centaurs), however, many H g ^> 9 objects are known from wide-field surveys, mostly 
detected at d < 20 AU. Thus the SO size distribution cannot remain at Of < 0.3. In 
fact, measurements of Jupiter Family Comet s (JFC s) in the H g £314-17 range give slopes 



Of 



0.5 ± 0.1 (see Table 6 of 



Solontoi et al. 



(120121 )). A divot can explain both a relative 



lack of objects beyond the break and the eventual recovery necessary to provide the JFCs. 
A divot also motivates the negative slopes measured, as a realistic divot will take the form 
of a decrease at the break, rather than the sharp discontinuity we use. We prefer the divot 
solution with ctf =0.5 and c ~6 (green star Fig. [6]) which matches the observations and 
allows for a single slope atf = 0.5 from the divot out to the H„ >14 Jupiter Fam ily comets 



whose slope is near the collisional equilibrium value (jO'Brien fc Greenberg 



20051 ). 
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5. External Arguments 



As the hot populations have similar colours and D > 100-km size distributions, 
i t seems likely th at they were all transplanted to join a pre-existing cold Kuiper Belt 



f lPetit et al 



201 ll ) during a common event early in the Solar System's history, and would 



thus logically share the same divot and small D distribution. Such a transpl ant process has 



been shown to successful 



2008 



Batvein et al. 



y implant TNOs in the stable Kuiper Belt regions ( iLevison et al. 



20111) . althou gh the resonant population ratios and i distribution are 



problematic ( IGladman et al 



20121 ). We thus look for evidence of such a feature in other 



hot populations. 



5.1. The Neptune Trojans 



A survey for Neptune Trojans ( ISheppard k, Truj illol 1 2 1 Oh provided significant evidence 
that ana ~ 0.8 power law cannot continue for D < 100-km Trojans; a divot was not 
apparent because Trojans significantl y smaller were not detected . The dispersed inclination 



distribution for the Neptune Trojans (ISheppard fc Trujillo 



20061), alt 



lough not yet precise ly 



measured, links thes e obje cts to all the other resonant populations (IGladman et al. 



20121 ). 



Sheppard fc Trujillo! (120101 ) used Neptune Trojan searches to argue that beyond ~23 
(corresponding to H g ~ 9) there was an absence of Trojans due to non-detections, and 
thus smaller Trojans were missing. Assuming that the Trojans and other resonant TNOs 
were implanted from a scattering population, and thus share the same size distribution, we 
confir med that our divot N(Hn) m atches the lack of D <100 km Neptune Torjan detections 



in the 



Sheppard fc Trujillo! (120101 ) surveys. 



Given our analysis, the conclusion would not be that small Neptune Trojans are 
"missing" , but rather that the sudden drop results in the population fainter than the divot 
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not recovering in on-sky surface density until at least H g > 11, by which point the deepest 
survey lacked the sensitivity to detect them. If correct, detection of several small (H g > 11) 
trojans would require surveying ~ 100 square degrees of sky to 26th magnitude at the 
correct elongation. 



5.2. hot Populations 



A recent deep telescopic survey (IFraser et al. 



20101 ) estimated a ~ 0.40 ± 0.15 (within 



error of our preferred «/ = 0.5) from the apparent-magnitude distribution for "close" 
(30 < d < 38) TNOs (orbits were not obtained). These distances are dominated by several 
hot populations, but the measurement is shallower tha n the u sual hot population slope of 



0.8. We calculated H„ magnitudes for the 



Fraser et al. 



(120 101 ) sample and find that due to 



the survey's great depth, this sample is dominated by H g > 9 TNOs and thus appears to 
measure the post-divot slope. 



6. Feasibility of a Divot 



A primordial size-distribution wave at small sizes (D = 2 km) could propagate (IFraser 



20091 ) to D ~100 km in a dynamically hot (Av = 2 km/s) collisi onal environment aft er 



500 Myr. Alternately , recent modeling of planetesimal creation (IJohansen et al. 



2007 



Morbidelli et al. 



20091 ). suggests that the protosolar nebula may only have produced 



planetesimals larger than a certain critical diameter, in which case the a& = 0.8 slope and 
the D ~ 100 km divot size are set by planetesimal formation physics; smaller objects 
appear only later due to collisional fragmentation. These scenarios match our results, 
where one interprets the hot population's N(H g ) to have been "frozen" when suddenly 
transplanted (scattered) from a denser region nearer the Sun to the large volume it now 
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occupies, ending the collisional evolution. An exciting prospect is that the divot directly 
records a preferential D that planet building produced in the solar-nebula region where 
the hot TNOs originally formed, and that the divot's depth (which could easily range from 
c=2-30) measures the integrated collisional evolution (depending on both the duration of 
the pre-scattering phase and th e random speeds present). An initi al distribution with no 



D <100 km TNOs was shown (ICampo Bagatin fc Benavidez 



20121 ) to evolve into a divot 



with c ~ 20 and ctf ~ 0.5, in the dy namical environment of the Nice model; such a 500-Myr 



quiescent phase ( IGomes et al 



20051 ) allows a divot to form but the evidence we find for 



a higher-i early phase may argue instead for a much shorter and more intense collisional 
environment. 

Our divoted N(H g ) produces a cumulative distribution (Fig. 2 B) with a shallow 



plateau for H„=9-12, similar to that deduced for the 



observations (IBernstein et al 



2004 : 



impactors of the saturnian moons ( IMinton et al 



Volk fc Malhotra 



lot p opulation and SOs in deep HST 



20081 ) and estimated for scattering 
20121 ). Our estimate of 2 x 10 6 SOs with 



H g < 13 and a slop e of at = 0.5 extrapolate s to ~ 2 x 10 9 SOs with H g <18, providing 



a sufficient number ( iDuncan fc Levison 



1997 



Volk fc Malhotra 



20081 ) of SOs to feed the 



Jupiter Family Comets, while satisfying the observed plateau. This novel parameterisation 
for the N(H g ) distribution (Fig. 2) matches our data set and simultaneously explains the 
puzzles of the JFC source, along with the "missing" Neptune Trojans, and the known 
rollover in the Kuiper Belt's luminosity function. 



Facilities: CFHT, NRC (HIA) 
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Fig. 1. — Cumulative distributions for 5 quantities, comparing the observed objects with 
the initially cold model's simulated detections. The black dash-dot curves show the orbital 
model's a, i, q, and d distributions; these differ from the detections (red staircase) due to 
observational biases. When coupled to two different N(H g ) distributions, the biases produce 
the differing predictions for the detections. A single power-law slope of a = 0.8 (blue dashed) 
is rejectable at > 99% confidence in d, q and H g . Contrastingly, our preferred divot N(H g ) 
(green curve, see Fig. 2) provides vastly better matches, although both produce too many 
\ow-i detections. 
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Fig. 2. — Histograms of the N(H g ) distribution for our preferred divot sol ution. The vertica l 



axes show the total SO numbers using the absolute CFEPS calibration ( jPetit et al. 



201 If ). 

A: The differential distribution (solid green), with an extrapolated a = 0.8 beyond H g = 9 
(dashed blue). The contrast c ~ 6 is the ratio of differential number on either side of the 
divot. B: The cumulative version. For H g >13 the cumulative N(< H g ) has reached slope 
a/ after the flattened region beyond the divot. 
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1.0 10.0 100.0 

contrast c 



Fig. 3. — Contour values (colors) computed from a grid (points) in (c, aj) space, giving 
the Anderson-Darling probability of drawing the most rejectable of the d, q, i, and H g 
distributions from the model of today's scattering objects, from a hot initial population (see 
text). White areas indicate when the worst distribution had <1% probability of coming from 
the model, and the black contour bounds <5%. A single power-law (blue star) is rejected. 
Our favoured model (green star) satisfies both (a) «/ ~ 0.5 like known JFCs, (b) a < 0.6 
which prevents the extrapolated mass of small SOs from diverging. 
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orbital inclination i (degrees) 



Fig. 4. — The hot (black dashed) and cold (light blue dashed) KRQ11 intrinsic i distributions 
produce the biased distributions (solid lines) for our preferred model for comparison with 
the CFEPS sample (red). The hot model significantly improves the match because it has 
a relative lack (at the current epoch) of low-z SOs to be detected by the survey; the eleven 
CFEPS SOs have < 1 % probability of being drawn from the initially cold simulation because 
it predicts i < 7° survey detections should form about half the detected sample (solid light 
blue curve). 
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Designation 


a (AU) 


e 


i (deg) 


d (AU) 


H, 


L4k09 


30.19 


0.18517 


13.586 


26.63 


9.5 


HL8al 


32.38 


0.37406 


42.827 


44.52 


7.3 


L4m01 


33.48 


0.33308 


8.205 


31.36 


8.9 


L4p07 


39.95 


0.28088 


23.545 


29.59 


7.7 


L3q01 


50.99 


0.484 


6.922 


38.17 


8.1 


L7a03 


59.61 


0.43949 


4.575 


46.99 


7.1 


L4vll 


60.04 


0.62928 


11.972 


26.76 


10.0 


L4v04 


64.10 


0.50638 


13.642 


31.85 


9.1 


L4vl5 


68.68 


0.6999 


14.033 


22.95 


9.0 


L3h08 


159.6 


0.7613 


15.499 


38.45 


8.0 


HL7j2 


133.25 


0.72374 


34.195 


37.38 


8.4 



Table 1: CFEPS actively scattering sample 



